Determination of the Newtonian Gravitational Constant Using Atom Interferometry 
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We present a new measurement of the Newtonian gravitational constant G based on cold atom 
interferometry. Freely falling samples of laser-cooled rubidium atoms are used in a gravity gra- 
diometer to probe the field generated by nearby source masses. In addition to its potential sen- 
sitivity, this method is intriguing as gravity is explored by a quantum system. We report a value 
of G = 6.667 • 10"^^ m'^ kg~^ s~^, estimating a statistical uncertainty of ±0.011 ■ 10^^^ m'^ kg^^ 
and a systematic uncertainty of ±0.003 • 10~^^ m'^kg"^ s~'^. The long-term stability of the instru- 
ment and the signal-to-noise ratio demonstrated here open interesting perspectives for pushing the 
measurement accuracy below the 100 ppm level. 
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The Newtonian constant of gravity G is one of the 
most measured fundamental physical constants and at 
the same time the least precisely known. Improving the 
knowledge of G has not only a pure metrological inter- 
est, but is also important for the key role that it plays 
in theories of gravitation, cosmology, particle physics, 
in geophysical models, and astrophysical observations. 
However, the extreme weakness of the gravitational force 
and the impossibility of shielding the effects of gravity 
make it difhcult to measure G, while keeping systematic 
effects well under control. Many of the measurements 
performed to date are based on the traditional torsion 
pendulum method Jj] , direct derivation of the historical 
experiment performed by Cavendish in 1798. Recently, 
many groups have set up new experiments based on dif- 
ferent concepts and with completely different systemat- 
ics: a beam-balance system [2|, a laser interferometry 
measurement of the acceleration of a freely falling test 
mass Q , experiments based on Fabry-Perot or microwave 
cavities However, the most precise measurements 

available today still show substantial discrepancies, lim- 
iting the accuracy of the 2006 CODATA recommended 
value for G to 1 part in 10*. From this point of view, 
the realization of conceptually different experiments can 
help to identify still hidden systematic effects and there- 
fore improve the confidence in the final result. 

Cold-atom interferometry has demonstrated outstand- 
ing performances for the measurement of tiny rotations 
and accelerations and it is widely used for many applica- 
tions: precision measurements of gravity Q, gravity gra- 
dient Q , and rotation of the Earth ^, 0] , but also tests of 
the Einstein's weak equivalence principle [l3| , tests of the 
Newton's law at short distances I 111, and measurement 
of fundamental physical constants |l2l. ITsj . Applications 
of these techniques for fundamental physics experiments 
in space are under study 

In this paper, we present a new determination of the 
Newtonian constant of gravity based on cold-atom in- 



terferometry. An atomic gravity gradiometer is used to 
measure the differential acceleration experienced by two 
freely falling samples of laser-cooled rubidium atoms un- 
der the influence of nearby tungsten masses. The mea- 
surement is repeated in two different configurations of 
the source masses and modeled by a numerical simula- 
tion. From the evolution of the atomic wavepackets and 
the distribution of the source masses, we evaluate the 
expected differential acceleration, having G as unique 
free parameter. A value for the Newton's constant of 
gravity is finally extracted by comparing experimental 
data and numerical simulations. Proof-of-principle ex- 
periments with similar schemes using lead masses were 
already presented in [ll,[3]. In the present work, specific 
efforts have been devoted to the control of systematic ef- 
fects related to atomic trajectories, positioning of source 
masses, and stray fields. In particular, the high density 
of tungsten is crucial in our experiment to compensate 
for the Earth gravity gradient. In both configurations of 
the source masses, atom interferometers can therefore be 
operated in spatial regions where the overall acceleration 
is slowly varying and the sensitivity of the measurement 
to the initial position and velocity of the atoms strongly 
reduced. 

Our atom interferometer uses light pulses to stimu- 
late *^Rb atoms on the two-photon Raman transition 
between the hyperfine levels F — 1 and F = 2 of the 



ground state [17|. The light field is generated by two 



counter-propagating laser beams with wave vectors ki 
and k2 ~ — ki aligned along the vertical axis. Atoms 
interact with the Raman lasers on the three-pulse se- 
quence 7r/2 — TT — 7r/2 which splits, redirects, and re- 
combines the atomic wavepackets. At the end of the 
interferometer, the probability of detecting the atoms 
in the state = 2 is given by P2 — 1/2 • (1 - cos$), 
where 4> represents the phase difference accumulated by 
the wavepackets along the two interferometer arms. In 
the presence of a gravity field, atoms experience a phase 
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FIG. 1: Schematic of the experiment showing the gravity gra- 
diometer set-up with the Raman beams propagating along the 
vertical direction. During the G measurement, the position 
of the source masses is alternated between configuration Ci 
(left) and C2 (right). 



shift $ = (ki — k2) • g depending on the local grav- 
itational acceleration g The gravity gradiometer 
consists of two absolute accelerometers operated in dif- 
ferential mode. Two spatially separated atomic clouds 
in free fall along the same vertical axis are simultane- 
ously interrogated by the same Raman beams to provide 
a measurement of the differential acceleration induced by 
gravity on the two samples. 

Figure [1] shows a schematic of the experiment. The 
gravity gradiometer set-up and the configurations of the 
source masses (Ci and C2) used for the G measurement 
are visible. Further details on the atom interferometer 
ap par atus and the source masses assembly can be found 
in 1^ 1^. A magneto-optical trap (MOT) with beams 
oriented in a 1-1-1 configuration collects rubidium atoms 
from the background vapor of the chamber at the bottom 
of the apparatus. Using the moving molasses technique, 
the sample is launched vertically along the symmetry axis 
of the vacuum tube and cooled down to a temperature 
of 2.5 ^K. The gravity gradient is probed by two atomic 
clouds moving in free fiight along the vertical axis of the 
apparatus and simultaneously reaching the apogees of 
their ballistic trajectories at 60 cm and 90 cm above the 
MOT. Such a geometry, requiring the preparation and 
the launch of two samples with high atom numbers in 
a time interval of about 100 ms, is achieved by juggling 
the atoms loaded in the MOT: A first cloud of 10® atoms 
is launched to 80 cm above the MOT; during its ballis- 
tic flight, 5-10® atoms are loaded and launched to 90 
cm; finally, the first cloud is recaptured in the MOT and 
launched again to a height of 60 cm. Shortly after this se- 
quence, the two atomic samples are velocity selected and 
prepared in the {F = l^mp = 0) level. The interferom- 



eters take place at the center of the vertical tube shown 
in Fig. [TJ In this region, surrounded by a system of two 
/z- metal shields (76 dB attenuation factor in the axial di- 
rection) , a uniform magnetic field of 250 mG along the 
vertical direction defines the quantization axis. The field 
gradient along this axis is lower than 5 • 10^^ G /mm. The 
population of the ground state is measured in the cham- 
ber placed just above the MOT by selectively exciting 
the atoms on the = 1, 2 hyperfine levels and detect- 
ing the light-induced fluorescence emission. We typically 
detect 10^ atoms on each rubidium sample after the in- 
terferometer sequence. 

The Raman beams are generated by two extended- 
cavity diode lasers and amplified in a tapered amplifier. 
An optical phase-locked loop keeps their frequency dif- 
ference in resonance with the transition between the two 
hyperfine levels of the ^^Rb ground state 6.8 GHz) and 
stabilizes their relative phase to about 100 mrad (rms) in 
a 5 Hz— 10 MHz bandwidth [i^]. The two co-propagating 
beams enter the vacuum system from the bottom of the 
MOT cell with the same linear polarization, travel along 
the axis of the vertical tube and, after crossing a quarter- 
wave plate, are refiected back producing a linJ-lin po- 
larization scheme. Their frequency difference is ramped 
down with continuous phase during the experiment cy- 
cle to compensate for the Doppler effect and ensure that 
freely falling atoms remain resonant with the laser light. 
In this way, only one pair of counterpropagating laser 
beams with crossed linear polarizations is able to stimu- 
late the atoms on the two-photon transition. The three- 
pulse interferometer has a duration of 2T — 320 ms. The 
TT pulse lasts 48 /xs and occurs 5 ms after the atomic clouds 
reach their apogees. Therefore, during the 7r/2 — tt — 7r/2 
pulse sequence, atomic trajectories are almost symmetric 
with respect to the apogees, allowing a better control of 
systematic effects induced by spurious magnetic fields. 

The source masses 1^ are composed of 24 tungsten 
alloy (INERMET IT180) cyhnders, for a total mass of 
about 516 kg. They are positioned on two titanium plat- 
forms and distributed in hexagonal symmetry around the 
vertical axis of the tube (see Fig. [T]). Each cylinder is 
machined to a diameter of 100.00 mm and a height of 
150.20 mm after a hot isostatic pressing treatment ap- 
plied to compress the material and reduce density inho- 
mogeneities. The internal structure of the material has 
been analyzed using different methods: ultrasonic tests, 
microscope analysis, and surface studies. Finally, a de- 
structive test has been performed to characterize the den- 
sity distribution. We have measured a maximum density 
variation of 2.6 • 10"'^ and an average density inhomo- 
geneity of 6.6 • 10~* over volume samples of about 1/40 
of the whole cylinder volume. 

To evaluate the stability of the gravity gradiometer, we 
performed a 5-hour measurement run in which the dif- 
ferential phase shift was continuously monitored. Each 
atom interferometer in the gravity gradiometer measures 
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FIG. 2: Acceleration along the symmetry axis of the system 
in configuration Ci (left) and C2 (right). The simulation in- 
cludes the efi'ect of the source masses and the Earth gravity 
gradient. The spatial extension of the atomic trajectories for 
the upper and lower interferometers is indicated by the thick 
lines. The shaded regions show the positions of the source 
masses. 



the local acceleration with respect to the common ref- 
erence frame identified by the wavefronts of the Raman 
lasers. Therefore, even if the phase noise induced by vi- 
brations on the retro-reflecting mirror completely washes 
out the atom interference fringes, the signals simulta- 
neously detected on the upper and lower accelerometer 
remain coupled and preserve a fixed phase relation. As a 
consequence, when the trace of the upper accelerometer 
is plotted as a function of the lower one, experimental 
points distribute along an ellipse. The differential phase 
shift is then obtained from the eccentricity and the ro- 
tation angle of the ellipse fitting the data [21| . We have 
evaluated the Allan variance of the differential phase shift 
and verified that it decays as the inverse of the integration 
time, showing the typical behavior expected for white 
noise. The instrument has a sensitivity of 140 mrad at 
1 s of integration time, corresponding to a sensitivity to 
differential accelerations of 3. 5 - 10^^ g in Is. These per- 
formances are compatible with a numerical simulation as- 
suming quantum-projection-noise-limited detection with 
10^ atoms. 

The Newtonian gravitational constant has been ob- 
tained from a series of gravity gradient measurements 
performed by periodically changing the vertical position 
of the source masses between configuration Ci and C2, 
with the atoms always launched along the same trajec- 
tories. Figure shows the acceleration along the sym- 
metry axis for both configurations. Because of the high 
density of tungsten, the gravitational field produced by 
the source masses is able to compensate for the Earth 
gravity gradient. Operating the interferometers close to 
these stationary points reduces the uncertainty on G due 
to the knowledge of the atomic positions by two orders 
of magnitude. A typical data sequence used for the mea- 
surement of G is shown in Fig. [3) Each phase measure- 
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FIG. 3: Typical data set showing the modulation of the differ- 
ential phase shift measured by the atomic gravity gradiometer 
when the distribution of the source masses is alternated be- 
tween configuration Ci (red squares) and C2 (blue squares). 
Each phase measurement is obtained by fitting a 24-point scan 
of the atom interference fringes to an ellipse. 



ment is obtained by fitting a 24-point scan of the atom 
interference fringes to an ellipse. The modulation on the 
differential phase shift measured by the gravity gradiome- 
ter can be resolved with a signal-to-noise ratio of 180 
after about one hour. Experimental data have been col- 
lected during two separate measurement runs performed 
in February and May 2007. 

The differential nature of the measurement ensures 
an efficient rejection of many systematic errors: Firstly, 
each gravity gradient measurement is highly insensitive 
to phase shifts appearing in common mode on the two 
interferometers [7|; in addition, the differential measure- 
ment performed by alternating the source masses be- 
tween configuration Ci and C2 removes effects not de- 
pending on the masses position. A budget of the sys- 
tematic errors affecting the value of G is presented in 
Table [H Each effect is taken into account in our nu- 
merical simulation and its contribution evaluated on the 
basis of the measurements performed on the relevant pa- 
rameters. Positioning errors account for uncertainties in 
the position of the 24 tungsten cylinders along the radial 
and vertical direction, both in configuration Ci and C2. 
Density inhomogeneities of the source masses have been 
modeled. An upper bound of the systematic uncertainty 
on G has been evaluated on the basis of the maximum 
density variation measured during the characterization 
of the tungsten cylinders. The contribution of the plat- 
forms holding the source masses has been included as 
well. The velocity of the atomic clouds and their po- 
sition at the time of the first interferometer pulse have 
been calibrated performing time-of-flight measurements 
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Systematic effect 


AG/G (xlO-*) 


Radiaf position 


1.2 


Vertical position in C'l 


2.7 


Vertical position in C2 


2.1 


Cylinders mass 


0.9 


Cylinders density inliomogeneity 


0.21 


Support platforms mass 


0.8 


Initial position of the atomic clouds 


0.18 


Initial velocity of the atomic clouds 


2.3 


Gravity gradient 


1 


Stability of the on-axis B-field 


0.3 


Stability of the launch direction 


0.6 


Total 


4.6 



TABLE I: Error budget of the systematic shifts affecting the 
G measurement. 



and detecting the atoms when crossing a horizontal light 
sheet twice, once on the way up and once on the way 
down. The knowledge of the gravity gradient is also im- 
portant to identify the best positions for the apogees of 
the two atomic clouds with respect to the source masses 
(see Fig. [2]). Finally, special care has to be taken of all the 
systematic errors that could depend on the source masses 
distribution. For instance, the presence of an inhomoge- 
ncous magnetic field along the vertical axis depending on 
the position of the source masses can introduce a system- 
atic phase shift affecting the G measurement. However, 
the low magnetic susceptibility of the tungsten alloy and 
the two /x-metal shields surrounding the interferometer 
region ensure an excellent control of stray fields and re- 
duce possible systematic errors below the measurement 
sensitivity of our apparatus. Using a tiltmeter, we have 
verified that the direction of the Raman lasers, which 
defines the sensitive axis of our interferometer, is stable 
within 1 /irad and does not depend on the position of the 
source masses. At the same time, the stability of the 
launch direction has been evaluated and an upper bound 
on possible phase shifts induced by the Coriolis acceler- 
ation on the differential measurement estimated. 

After an analysis of the error sources affecting our 
measurement, we obtain a value of G = 6.667 ■ 
10""'^^ m'^ kg~^ s~^, with a statistical uncertainty of 
±0.011 • 10"^"'^ m^ kg""'^ s~^ and a systematic uncertainty 
of ±0.003 ■ 10"" m^ kg"^ s-2_ FigureHshows this result 
and compares it to the most precise G values recently 
obtained in different experiments. Our measurement is 
consistent with the 2006 CODATA value within one stan- 
dard deviation. 

We have reported a new measurement of the Newto- 
nian constant of gravity based on atom interferometry. 
This technique, completely different from the methods 
used to date, explores new systematics, globally improv- 
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FIG. 4: Our result of the Newtonian gravitational constant 
compared to the most precise G measurements recently ob- 
tained and to CODATA recommended values. 



ing the confidence in the knowledge of G. Presently, our 
statistical error averages down reaching a value compa- 
rable to the systematic uncertainty in about one day of 
integration time. As expected, the main contribution to 
the systematic error on the G measurement derives from 
the positioning accuracy of the source masses. This er- 
ror will be reduced by about one order of magnitude once 
the position of the tungsten cylinders will be measured 
by a laser tracker. From the analysis of systematic effects 
and the characterization of the stability of our appara- 
tus, we foresee interesting perspectives for pushing the 
measurement accuracy below the 100 ppm level. 
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